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Abstract 

We show a numerical evidence that a tethered surface model with extrinsic curvature 
undergoes a first-order crumpling transition between the smooth phase and a non- 
smooth phase on triangulated tori. The results obtained in this Letter together with 
the previous ones on spherical surfaces lead us to conclude that the tethered surface 
model undergoes a first-order transition on compact surfaces. 
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1 Introduction 



The tethered surface model of Helfrich, Polyakov, and Kleinert (HPK) [1,2,3] 
is a surface model for artificial and biological membranes and for strings 
[4,5,6,7,8]. The vertices, which represent the lipid molecules of membranes, 
are prohibited from the free diffusion over the surface because of the fixed 
connectivity nature of the triangulated surface. The phase transition between 
the smooth phase and a non-smooth phase (or crumpled phase) of the model 
has long been studied, and the phase structure is being clarified step by step 
[14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32]. 

An intrinsic curvature model for membranes is known to undergo a first-order 
crumpling transition, which was confirmed by Monte Carlo (MC) simulations 
[15,16,17]. The transition occurs on a disk [16], which is a noncompact surface. 
Moreover, the transition can be seen on a sphere [15] as well as on a torus 
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[17]. Therefore, we can state that the tethered surface model with intrinsic 
curvature has a first-order transition on an arbitrary surface in R 3 . 

However, it remains to be stated that the extrinsic curvature model has a first- 
order transition on compact surfaces. In fact, the first-order transition of the 
model has never been confirmed on a torus, although it has been confirmed 
on a sphere by the MC simulations [25,26,27]. The reason why it is so hard to 
clarify the order of the transition on a torus seems because of the size effect. 
The extrinsic curvature model shows the first-order transition only on the 
sphere of size iV>7000 [26,27], which is very large in contract to that in the 
model with intrinsic curvature, where the transition can be seen even on very 
small surfaces [15,16,17]. 

It is true that the size effect disappears at sufficiently small curvatures on the 
sphere because the first-order transition occurs only at large N, which corre- 
sponds to the large radius R of the sphere. Because the principal curvature is 
given by k=1/R, we have k = \ where the area of the surface is identi- 
fied with N. On the torus in Fig. 1(a) made from the square of size (Li,!^) 

shown in Fig. 1(b) we have approximately k = ■J^jf, which is the larger one 
of the two principal curvatures, where L\ = 2L2 and lisRi = Li(i = 1, 2) were 
assumed. Therefore, the first-order transition can be seen on the torus of size 
6 ~ 7 times larger than that of the sphere, because the sphere and the torus 
are considered to have the common k where the size effect disappears. It is 
natural to consider that low-frequency modes of surface fluctuations play an 
important role in the crumpling transition, and that the corresponding energy 
scales are independent of the surface topology. 




In this Letter, we study the extrinsic curvature model in [25], where the first- 
order transition was seen even on the sphere of size N = 1500 ~ 2500, which are 
relatively smaller than iV ~ 7000 where the first-order transition was observed 
in the model with the standard bending energy. The bending energy in [25] 
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is slightly different from the standard one, however, the model is in the same 
class as that of the standard bending energy. 

We will show in this Letter the first numerical evidence that a tethered surface 
model with extrinsic curvature undergoes a first-order crumpling transition 
between the smooth phase and a non-smooth phase on triangulated surfaces 
of torus topology. It was already reported by us that the tethered surface model 
of HPK undergoes a first-order transition on spherical surfaces [25,26,27]. The 
results obtained in this Letter together with the previous ones [25,26,27] lead 
us to conclude that the tethered surface model of HPK undergoes a first-order 
crumpling transition on compact surfaces. 



2 Model 



In order to define the bending energy S 2 , we use the normal vector of the 
vertex i such as 



N- 



where J2j(i) denotes the summation over triangles linked to the vertex i. 
The vector is the unit normal of the triangle and Aa j{i) is the area 
of 

The Gaussian tethering potential Si and the extrinsic curvature S 2 are defined 
by 

S 1 = '£(X i -X j ) 2 , S 2 = ££[l-n(0-n jW ], (2) 

(ij) » 



where Yluj) is the sum over bond (ij) connecting the vertices and Xj. It 
should be noted that S^illdef) = (1 ~~ n W ' n (j)) defined only by using 
the normal vectors n(i) in Eq.(l) is not well defined. This ill-definedness comes 
from the fact that there exist two surfaces locally different from each other 
that has the same value of S^illdef). Two normal vectors at the ends of a 
bond can be parallel for surfaces that are not smooth. 

In order to understand graphically the interaction described in Eq. (2), we 
show the interaction range between the normals of triangles in Fig. 2. The 
normal of the shaded triangle interacts with the normals of the triangles in 
the figure. 

The partition function is defined by 
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Fig. 2. The range of the interaction between the normal vectors in S2. The normal of 
the shaded triangle interacts with the normals of the triangles shown in the figure. 



where N is the total number of vertices as described above. The expression 
S(X) shows that S explicitly depends on the variable X. The coefficient b is 
the bending rigidity. The surfaces are allowed to self-intersect. The center of 
surface is fixed in Z(b) to remove the translational zero-mode. 



3 Monte Carlo technique 

The canonical Monte Carlo technique is used to update the variables X so 
that X' = X + 5X , where the small change SX is made at random in a small 
sphere in R 3 . The radius 5r of the small sphere is chosen at the beginning 
of the simulations to maintain the rate of acceptance rx for the X-update as 
r x — 0.5. 

We use a random number called Mersenne Twister [33] in the MC simulations. 
Single sequence of random number is used for 3-dimensional move of vertices 
X and for the Metropolis accept/reject in the update of X. 

We use surfaces of size N = 1800, N = 3200, N = 5000, N = 9800 and N = 16200 
for Li/L 2 = 2, and N = 1764, N = 3600, N = 6400, N = 10000 and N = 16900 
for L1/L2 = 4. The size of the surfaces is obviously larger than that of the 
same model on spherical surfaces used in [25]. 

The total number of MCS after the thermalization MCS is about 2 x 10 8 ~ 
2.6 x 10 8 at the transition point of surfaces of size iV = 9800 ~ N— 16900, and 
1.0 x 10 8 ~ 1.6 x 10 8 for N = 1764 ~ N = 6400. Relatively smaller number of 
MCS is done at non-transition points in each N. 




i=l 



(3) 



S(X) - Si + 65 : 



2- 
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4 Results 
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Fig. 3. Si/N vs. b obtained on surfaces of (a) JV = 5000 ~ 16200, L 1 /L 2 = 2, and 
(b) iV = 6400 ~ 16900, L 1 /L 2 = A. 

The relation S\/N =1.5 is expected from the scale invariance of the partition 
function and can be used to check that the simulation was correctly performed. 
We also expect that this relation should not be influenced by whether the phase 
transition is of first order or not. Figures 3(a) and 3(b) are Si/N against b 
obtained on surfaces of size N = 5000 ~ 16200, Li/L 2 = 2 and N = 6400 ~ 
16900, Li/L 2 = 4. We find from these figures that the expected relation S±/N = 
1.5 is satisfied. 
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Fig. 4. X 2 vs. b obtained on surfaces of (a) iV = 5000 ~ 16200, L 1 /L 2 = 2, and (b) 
N = 6400 ~ 16900, L 1 /L 2 = A. 

The surfaces become smooth at b^oo and crumpled at b^O. Hence the size 
of the surfaces can be reflected in the mean square size X 2 , which is defined 
by 



(4) 



In order to see the size of surfaces, we plot X 2 against b in Fig. 4(a). We 
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see that X 2 grows rapidly at intermediate b, however it seems continuously 
change against b. We must note that the continuous behavior of X 2 does not 
always indicate a higher-order transition between the smooth phase and the 
crumpled phase. It is possible that X 2 changes continuously at the transition 
point, because the phase transition is characterized by the analytic structure 
of the partition function: the transition is of first-order (second-order) when 
|j| is discontinuous is discontinuous). 




0.455 u 0.465 0.455 K 0.465 

(a) b (b) b 

Fig. 5. The bending energy S2/2NB vs. b obtained on surfaces of (a) 
N = 5000 ~ 16200, Li/L 2 = 2, and (b) N = 6400 ~ 16900, Li/L 2 = 4. N B is 
the total number of bonds. 

The bending energy S2/2NB is plotted in Figs. 5(a),(b). By dividing S2 by 
2Nb, we obtain the bending energy per bond. In fact, the definition of 62 in 
Eq. (2) includes the summation over bonds twice. S2/2NB rapidly changes at 
b where X 2 rapidly changes, however, discontinuous change is also hard to see 
in S 2 /2N B . 
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Fig. 6. The specific heat Cs 2 vs. b obtained on surfaces of (a) N = 5000 ~ 16200, 
Li/L 2 = 2, and (b) iV = 6400 ~ 16900, Li/L 2 = 4. 

The specific heat C S2 is defined by 

C, 2 = ^((S 2 -(S 2 » 2 >, (5) 
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which reflects the phase transition if it anomalously behaves. Figure 6(a),(b) 
shows Cs 2 against b, which were obtained on surfaces of N = 5000 ~ 16200, 
L 1 /L 2 = 2, and iV = 6400 ~ 16900, L 1 /L 2 = 4, respectively. Anomalous behav- 
iors are obviously seen in C S2 in the figures. This indicates that the fluctuation 
of S 2 becomes very large at b where Cs 2 has the peak at each N. Thus, we 
found a sign of a phase transition. 




Fig. 7. Log-log plots of the peak value C<? ax against N, which were obtained on 
surfaces of (a) iV = 1800 ~ 16200, L 1 /L 2 = 2, and (b) iV = 1764 ~ 16900, L 1 /L 2 = 4. 

In order to see an order of the transition, we plot the peak value C!? ax against 
N in Figs. 7(a), (b) in a log-log scale. Figures include the simulation results 
obtained on surfaces of TV = 1800, N = 3200 for L x jL 2 = 2, and on iV= 1764, 
N = 3600 for Li/L 2 = 4. The straight lines in Figs. 7(a), (b) were drawn by 
fitting the largest three C™ 2 ax to 

C™ ax ~ N a , (6) 



where a is a critical exponent of the transition. Thus, we have 



o-i = 1.065 ± 0.090, (Li/L 2 = 2), 

a 2 = 1.018 ±0.082, {L x jL 2 = {). (7) 

These results indicate that the phase transition is of first-order in both cases 
for Li/L 2 = 2 and L 1 /L 2 = 4. We also see that the model show a continuous 
transition on small surfaces of size iV < 5000 ~ 7000, since the exponent a is 
a < 1 at N < 5000 ~ 7000, which is obvious from the figures. Therefore, the 
torus is severly influenced by the size effect as expected, since the first-order 
transition can be seen in the same model on the sphere of N > 1500 ~ 2500 
[25]. 

Figures 8(a), (b) are the autocorrelation coefficient A(S 2 ) of £2 defined by 
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Fig. 8. Auto-correlation coefficient A(S2) obtained on surfaces of (a) N = 16200, 
Li/L 2 = 2, (b) N = 16900, Li/L 2 = 4, and A(X 2 ) obtained on surfaces of (c) 
N = 16200, L 1 /L 2 = 2, (d) N = 16900, Li/L 2 = 4. 



A(5 2 ) = 



E,52(r,)5 2 (T j+1 ) 



Tj + i = Tj + n x 500, n — 1,2, 



(8) 



The autocorrelation coefficient A(X 2 ) of X 2 were also plotted in Figs. 8(c),(d). 
The phase transition can be reflected in the convergence speed of the MC 
simulations. The horizontal axes in the figure represent 500 xn (n — 1,2, • • •)- 
MCS, which is a sampling-sweep between the samples S 2 (Ti) and S 2 (T i+i ) . 
The coefficients A(5' 2 ) and A(X 2 ) in Figs. 8(a)-(d) were obtained on surfaces 
of size N = 16200 for L x jL 2 = 2 and N = 16900 for L x jL 2 = 4. The symbol 
(O) represents data obtained at the transition point, and the symbols (A) 
and (y) represent those obtained in the crumpled phase and in the smooth 
phase, respectively. 

We clearly see from the figures that the convergence speed at the transition 
point (O) is about 10 times larger than those in the crumpled (A) and the 
smooth (v) phases both for S 2 and X 2 . This is a phenomenon of the critical 
slowing down typical to the phase transition. The reason why the convergence 
speed at the transition point is so larger than in the crumpled and the smooth 
phases is that the phase space volume (C R 3 ), where the vertices X, take 
their values, becomes larger at the transition point than that in the smooth 
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and crumpled phases. 



(a) Smooth surface (b) The surface section of (a) 

iV = 16200, Li/L 2 = 2, 
6 = 0.456 




(c) Crumpled surface (d) The surface section of (c) 

iV = 16200, Li/L 2 = 2, 
6 = 0.456 

Fig. 9. Snapshot of the surface of iV = 16200, L\/L2 = 2 obtained at the transition 
point b = 0.456; (a) a smooth surface and (b) the surface section, (c) a crumpled 
surface and (d) the surface section. The mean square size of the smooth surface is 
X 2 = 87, and that of the crumpled one is X 2 = 32. 

Snapshots of surfaces are shown in Figs. 9(a)-(d), which were obtained at the 
transition point 6 = 0.456 on the N = 16200, L\jL2 = 1 surface. The surface in 
Fig. 9(a) is a smooth one, whose X 2 is about X 2 = 87. Figure 9(b) is the section 
of the surface in Fig. 9(a). The surface in Fig. 9(c) is a crumpled one, whose 
X 2 is about X 2 = 32. Figure 9(d) is the section of the surface in Fig. 9(c). We 
can see from the surface section in Fig. 9(b) that there is a lot of empty space 
inside the smooth surface in Fig. 9(a). On the contrary, we can also see from 
Fig. 9(d) that the space is obviously reduced in the crumpled surface in Fig. 
9(c). We note that X 2 seems continuously change at the transition point, in 
contrast to the model on a sphere [25,26,27], where X 2 has two distinct values 
at the transition point. The variation of X 2 against MCS, which was not 
depicted in a figure, indicates the continuous behavior of X 2 at the transition 
point. 
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(a) Smooth surface (b) The surface section of (a) 

iV = 16900, Li/L 2 = 4, 
6 = 0.456 




(c) Crumpled surface (d) The surface section of (c) 

iV = 16900, Li/L 2 = 4, 
6 = 0.456 

Fig. 10. Snapshot of the surface of N = 16900, Li/L2 = 4 obtained at the transition 
point b = 0.456; (a) a smooth surface and (b) the surface section, (c) a crumpled 
surface and (d) the surface section. The mean square size of the smooth surface is 
X 2 = 98, and that of the crumpled one is X 2 = 31. 



Snapshots of surfaces are shown in Figs. 10(a)— (d), which were obtained at 
the transition point 6 = 0.456 on the N = 16900, Li/L2 = 4 surface. The surface 
in Fig. 10(a) is a smooth one, whose X 2 is about X 2 = 98. Figure 10(b) is the 
section of the surface in Fig. 10(a). The surface in Fig. 10(c) is a crumpled one, 
whose X 2 is about X 2 = 31. Figure 10(d) is the section of the surface in Fig. 
10(c). We can see from the surface section in Fig. 10(b) that there is a lot of 
empty space inside the smooth surface in Fig. 10(a). On the contrary, we can 
also see from Fig. 10(d) that the space is obviously reduced in the crumpled 
surface in Fig. 10(c). We note also that X 2 seems continuously change at the 
transition point, as in the case for L 1 /L 2 = 2. 
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5 Summary and conclusion 



We have shown that a tethered surface model with an extrinsic curvature 
undergoes a first-order phase transition between the smooth phase and the 
crumpled phase on tori. The tori were constructed by connecting the sides 
of rectangular surfaces of size L\ x L 2 . The specific heat C$ 2 has the peak 
C!? ax at the transition point, and C!? ax scales according to C!? ax ~ N a , where 
a = 1.065(90) for Li/L 2 = 2 and a = 1.018(82) for L x jL 2 = 2. These results 
indicate the transition is of first-order on the basis of the finite-size scaling. 
We found an anomalous behavior of Cg^ x , which is consistent with the first- 
order transition on the surfaces of size N > 5000 ~ 7000. The surface size 
N = 5000 ~ 7000 is very large compared with N = 1500 ~ 2500 of the same 
model on a sphere [25]. This implies that very large tori are necessary for the 
model with the standard bending energy to undergo the first-order transition, 
because the model with the standard bending energy on a sphere shows the 
first-order transition on the surfaces of N> 7000 [26,27]. 

It must be noted that a gap was not seen in S 2 at the transition point, and 
hence the histogram of S 2 has no double peak structure in contrast to the 
model on spheres. The reason of this seems due to the size effect. The model 
on the torus is severely affected by the size effect. Nevertheless, we confirmed 
from the finite-size scaling that the model undergoes a fist-order transition as 
stated above. 

Thus, the results of the extrinsic curvature model on tori in this Letter together 
with those previously obtained on the same class of model on spheres can 
lead us to conclude that the model undergoes a first-order transition between 
the smooth phase and the crumpled phase on compact surfaces. Moreover, 
the phase transition and the topology change is compatible in the extrinsic 
curvature model. The term compatible means that both of two phenomena lead 
to the same result without depending on which phenomenon firstly occurs, as 
discussed in [17] for the intrinsic curvature model. 

It was already reported that the intrinsic curvature model undergoes a first- 
order transition, which is independent of whether the surface is compact or 
noncompact. This fact and the result in this Letter imply that the tethered 
surface model has a first-order crumpling transition on compact surface, and 
the transition occurs independent of whether the curvature energy is an in- 
trinsic one or an extrinsic one. 

This work is supported in part by a Grant-in- Aid for Scientific Research, No. 
15560160. 
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